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Abstract. We derive WKB approximations for a class of Airy and parabolic cylin- 
der functions in the complex plane, including quantitative error bounds. We prove 
that all zeros of the Airy function lie on a ray in the complex plane, and that the 
parabolic cylinder functions have no zeros. We analyze a limiting case where the 
parabolic cylinder functions go over to the Airy function. 
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1. Introduction and Statement of Results 

This paper represents an important first step in our program of proving the linearized 
stability of the Kerr back hole. The mathematical problem can be framed as showing 
that solutions of the so-called Teukolsky equation |1Q[ Q] for spin s=2, with smooth 
compactly supported initial data outside the black hole, decay uniformly on compact 
sets. The Teukolsky equation separates into an angular and radial ODE. The angular 
equation involves a Sturm-Liouville operator of the form 

Here the potential W is complex, and thus the operator is not self-adjoint. Following 
the procedure developed in [2], we seek a spectral representation for this angular op- 
erator. Our method for proving such a spectral representation requires detailed infor- 
mation on the eigensolutions. This information can be obtained by "glueing together" 
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approximate solutions and controlling the error using the methods in [3] . One method 
for obtaining such approximate solutions is to solve the Sturm-Liouville equation 

.^ + V{u))<t>[u) = (1.1) 

where V is a linear or quadratic polynomial with complex coefficients. The correspond- 
ing solutions are Airy and parabolic cylinder functions, respectively. In this paper we 
analyze properties of these special functions. 

Both of these special functions have interesting representations as contour integrals 
in the complex plane. Since it is difficult to estimate these integral representations 
directly, we analyze them with stationary phase-like methods (for an introduction to 
these methods see for example Section 7.7]) to obtain approximate WKB solutions 

J^expfTw) (1.2) 



with quantitative error bounds. In order to apply the methods in [3], we also need to 
control the function y := <j>'/(j>, which solves the associated Riccati equation 

y' = V-y 2 . (1.3) 

For y to be well-behaved, we must show that cf> has no zeros. This motivates our 
interest in ruling out zeros of the Airy and parabolic cylinder functions. 

Apart from being of independent interest, the results obtained here will be used in 
the forthcoming papers [HE]. 

2. Estimates for Airy Functions 

In this section, we assume that V is a linear function, 

V{u)=a + bu with a, b G C . (2.1) 

Then the Sturm-Liouville equation (|l.ip can be solved explicitly in terms of Airy 
functions, 

4>{u) = 2l(z) where z = b'i{a + bu) (2.2) 
and 21 is a linear combination of the Airy functions Ai and Bi (see [8j) 

2t(z) := -i^ (Ai(z) + zBi(z)) 

(the reason for our specific linear combination is that it has a particularly simple 
WKB asymptotic form; see (|2.5p below). The corresponding solution of the Riccati 
equation (II. 3p is given by 

= M 

Our first goal is to get asymptotic expansions and rigorous estimates of the func- 
tions 21' (z) and 2t(z) (see Theorem 12.11 below) . 

Using the integral representations in [8j eqns (9.5.4) and (9.5.5)] 

^ /-Doe™/ 3 / t 3 \ 

A K Z ) = 77— / exp( — - zt) dt 

Bi(*) = - J ^ exp ^- - zt) dt + - j ^ exp ^- - zt) dt 
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(where ooe* 7 ' refers to the end point of the contour lim t ^ +00 te lip , etc.), we obtain 

%{z) = -LJexp(j- zt\ dt , (2.4) 



where T is the contour T = — R + U e 7 ™ / ' 3 R + . Note that the last integral is finite 
because the factor e* / 3 decays exponentially at both ends of the contour. From the 
computation 

*3 



2l"(z) = -^= j t 2 exp - zt\ dt 



SX z+ ^)^ v ^~ zt ) dt=z%{z) 



1 /" / d \ /i 3 
vWr 

one immediately verifies that 21 is a solution of the Airy equation. 



2.1. WKB Estimates. We expect that for large z, the Airy solution should go over 
to the WKB wave function (jl.2p corresponding to the Airy potential V{z) = z, 



1 2 -A 



21wkb(*) :=z-*e* zi . (2.5) 

Here we define the roots by z a = exp(a log(z)), where the logarithm has a branch cut 
along the ray 

e"~ R+ . (2.6) 

In the next lemma, we show that with this branch cut, the WKB wave function ()2.5p 
approximates the Airy function 21 for large z, and we also derive rigorous error bounds. 

Theorem 2.1. Assume that for given e G (0, 5), 



Then 



2l(z) 



21 W kb(^) 

*'(*)- *Wb(*) 



argz0(-y-e, — + mod 2-7T . (2.7) 



< 7~T ■ J ^ (l + l°g 2 (4sin(e/2)|z|l) ) (2..S) 



z 



% sin J (e/2) 



21 W kb(^) 



^ — . o, 1 7TT (l + log 2 (4sin(e/2)|z|5 



2 



3 sin 2 (e/2) 



1 1 

|ip sin 3 (e/2) 



+ 7^2 ■ 3/./^ (l+|log(4sin(g/2)|z|l)| 3 ) . (2.9) 



Proof. The assumption (12. 7ft and our branch convention for the square root imply that 
the parameter c defined by 

c = Re(e"f (2.10) 
is positive and bounded by (see Figure [I]) 

c > \z\% sin(e/2) . (2.11) 
Rewriting the argument of the exponential in (|2.4p as 

i 3 2 3 r 3 

zt = - z? — ^fz t 2 H with r := t + \fz , 

3 3 3 
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we obtain 



-iiv/6 "| 



'2 3 
CX]) I - Z2 



e e 



exp ( - z2 



e 3 ^ z s 



(It 



(2.12) 
(2.13) 



where in f)2. 12|) we deformed the contour, whereas in the last line we introduced the 
new integration variable s = e J7T//6 r. Note that in (|2. 13[) the factor e~ ls / 3 is merely a 
phase, but the quadratic term in the exponent still ensures convergence of the integral 
in view of (pOOD . 

If we drop the cubic term in the exponent, the resulting Gaussian integral can be 
computed to give precisely the WKB wave function, 



21wkb(^) 



1 



7T 



'2 3 

exp | - z 2 — e 3 



' z s 



ds 



We thus obtain the error term 

E(z) :=2l(z)-2l W KB^) 



= — e « / 
V 71 " Jk 

After estimating the error term by 



exp [ - z 2 — e 



zs 2 ) e ~§ 



1 ds 



< / e — ' 


IS 3 

e 3—1 


ds 


J — oo 







(2.14) 



(2.15) 
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we decompose the integral into integrals over the regions [-L, L] and R\ [-L, L], which 
we bound as follows, 



-L,L] 



e 3 



1 



ds < 2 



ds 



< 2 e 2" 



-L,L] 
2 roo 



e 2 ds 



2^/27? _ £ i? 



e 2 



/: 



e _ — - 1 



L g 3 L A 

ds < 2 / — < — 
'o 3 " 6 



(2.16) 



(2.17) 



where in the last line we used that for all x G R the inequality |e* x — 1| < \x\ holds. In 
order to make both errors of about the same size, we choose 



L 



^ log 2 ((4 



3 . 
C)2 



(2.18) 



We thus obtain 



9 3 

— 2 2 2 

e s 



< 



2-§ 



1 4 



- + 67^ ? log2(4c)3/2 - ? ^ + log2(4c) ) • (2 - 19) 

Applying (l2~TT]l and (|33]| gives (pl8l) . 

In order to derive (|2.9j) , we differentiate (|2.13p and (|2.14p to obtain similar to (|2.15p , 
the identity 



(&'(*) -*Wb(*)) 



e e 



1 £7T rf\ ( 2 3 17T , p 

z — - e 3 s j exp l-z2— e -a z s 



Thus 



e-l^(j-z{K'{z)-^ KB (z))-zE(z)) <-i= [ 



s 2 e~ cs 



e 3 



e 3 



1 ds. 



ds 



Estimating the integral exactly as in (|2.16p and (|2.17p . we obtain 



2v^F 



s 2 e~ csZ 



e 3 



1 _ct£ 1 L 6 

ds < —5- e 2-1 _ — 

~ cf 20F 9 



Choosing L again according to (|2.18p . we apply (|2.19p to obtain 



„ 3 



^ (*'(*) -OWkb(*)) 



< 



e~§ J z£(z) +4.(l + |log(4c)| 3 ) 



1 



< O (1 + log 2 (4c)) + - (1 + I log(4c)| 3 ) 



C 2 1 ' " C 

The proof is completed by applying (|2.5|) and using (|2.11|) . 



□ 



In particular, this theorem allows us to take the limit as z goes to infinity along a 
ray through the origin, as long as we stay away from the branch cut (|2.6p . 

Corollary 2.2. Let z G C be off the ray e"fl+. Then 

_I %f(tZo) — 
t->oo 2l(t^o) 

where we again used the sign convention for the square root introduced just before (|2.6p . 
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Proof. We choose e so small that zq satisfies the condition (|2.7p . Then 

lim H ^ = lim H • 
and computing $i' WKB h° m t|2.5|) gives the result. □ 

2.2. An Estimate for the Riccati Equation for a Complex Potential. Our next 
goal is to show that the Airy function 21 has no zeros except on the branch cut (|2.6p . 
In preparation for this, we now derive an estimate for solutions of the Riccati equation 
for a general potential V with ImF > 0. Thus let y be a solution of the Riccati 
equation ([1.3|) . Decomposing y into its real and imaginary parts, y = a + i/3, we 
obtain the system 

a' = ReV-a 2 + p 2 (2.20) 
p' = lmV-2a(3. (2.21) 

Moreover, we set 

a = exp ^2 J cij . (2.22) 

Let us assume that ImV > on an interval [ttc^i] and that (a,/3) is a solution on 
a closed subinterval [^0,^2] C [uo,ui]. If /3(«o) > 0, it follows immediately from (|2.2ip 
that (3 > on the whole interval [uo,^]- A short calculation yields 

(o-/3)'(u) = almV>0. 

Hence the function a(3 is monotone increasing and 

(3(u) > /3 a —3— , 
a{u) 

where (3q = (3(uq) and similarly for all other functions. Using this inequality in (|2.20p . 
we obtain 

a'>ReV-a 2 + ^l, (2.23) 
a{u) z 

and integrating this inequality gives the following estimate. 

Lemma 2.3. Assume that ImV > on the interval [uq,ui] and that (a, j3) is a 
solution on a closed subinterval [uq,U2[ C [uq,u\] with P{uq) > 0. Then the function a 
satisfies the upper bound 

a(u) < g(u) on [u ,u 2 } , 

where we set 



and 



'C sinh(x) + a(uo) cosh(x 



cosh(x) + a{u2) sinh(x)/vC 
C= inf ReUGlR, x = VC(u-u 1 ). 

[Uo,Ul] 



Proof Dropping the last summand in (12.23[) . we obtain the inequality 

a' >C-a 2 . 

Then a < g, where g is the solution of the corresponding differential equation 

g' = C - g 2 , g(u 2 ) = a(u 2 ) . 
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Solving this differential equation gives the result. □ 
Using the result of this lemma in (|2.22p . we get 

a(u) = a exp I 2 / a) < a exp 2 / g J . 

\ Juq / \ Juq / 

Computing the integral, we obtain 

/-v/Ccosh iy/C{u - n 2 )) + a(u 2 ) sinh (y/C{u - u 2 )) \ 

a[u) < CTo — j= -^—j= r i—j= -r I . 

y V C cosh (v C(n — no)) + a(«2) sinh (vC(w- no)) / 

Using this estimate in (|2.23[) and setting u = u 2 gives the following result. 

Proposition 2.4. Assume that ImV > on [uq,u\] and that /3(uq) > 0. As- 
sume furthermore that the solution (a,/3) exists on the interval [uo,u] C [no,ni] and 
that a'(u) < 0. Then 



> C - a(n) 2 + ^ 



cosh (v / C(n — no)) H -1= sinh (v / C(n — no)) 

vC 

Corollary 2.5. Assume that either IvaV\< UQ>Ul ] > and (3{uq) > or ImV|[ U0)Ul ] < 
and /3(uq) < 0. T/ien the solution (a,/3) exists and is bounded on the interval [uq,ui\. 

Proof. It suffices to consider the case ImU|[ UOjU1 ] > and /3(uq) > because the other 
case is obtained by taking the complex conjugate of the Riccati equation (jl.3D . 

Let us assume conversely that the solution (a,/3) blows up at a point u 2 £ [no,ni]. 
We can assume that a blows up, because otherwise /3 could be obtained by integrat- 
ing $zm, 

P(u) = e L o a (3(u ) + / e- 2 I>ImV(T)dT . (2.24) 

J uo 

Let us assume that a tends to +oo. Then (|2.24|) shows that j3 stays bounded, and thus 
the right hand side of (I2.20p tends to — oo, a contradiction. On the other hand, if a 
is unbounded from below, then there is a sequence v n G [no,n2) with a(v n ) — > — oo 
and a'(v n ) < 0. This contradicts Proposition 12.41 □ 



2.3. Ruling out Zeros of the Airy Function. Our method for ruling out zeros of 
the Airy function 21 at a given point zq E C is to consider the Airy function along 
straight lines of the form 

z(u) = z + An with n € M , (2.25) 

where A is a complex parameter. Then the function (j>(u) := 2l(z(n)) satisfies the 
Sturm-Liouville equation (II. lj) with 

V(u) = A 2 (zq + An) . 

Applying Corollary 12.51 to the corresponding Riccati solution y := 4>'/<p yields the 
following result. 

Theorem 2.6. The function 21 has no zeros off the ray e~ ?7r//3 R + . 
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Figure 2. Admissible phases. 



A = e ltp with if e ( J, |) U (J, J) ; (2.26) 



Proof. Let zo be a point off the ray e l7r / 3 M + . We choose A such that (—A) is also off 
the ray e _l?r//3 M + . Then Corollary 12.21 applies and yields 

v I i— A / \ , r , ,-1 , A-r 

lim \u\ % y{u) = A lim u 2 zti = A v — A . 

u^-oo u— >— oo 2l(z(n)) 

Using our sign convention for the square root, we may choose A such that 

'IT 7T\ /7T 7T 

.6' 3 J U v3' 2 
then 

lim |u|-3 y( u ) = -e 3iip (2.27) 
it— y— oo 

(Note that the parameter —A = e 2%ip is indeed off the ray e _J7r//3 M + because tt/3 is 
excluded in (|2.26p .) Moreover, 

Im V(u) = Im(X 2 z ) + u Im(A 3 ) = Im(e 4 "%) - u Im(e 6 ^) . (2.28) 

In the considered range of if, the imaginary parts of e 3l<f and e 6lv have opposite signs 
(see Figure [2]). As a consequence, one sees from (|2.27p and (|2.28p that for sufficiently 
large negative u, the functions Imy and ImV are non-zero and have the same sign. 
Moreover, we can choose tp such that the imaginary parts of e ^zq and e GlLp have the 
same sign (see Figure [2]) . We conclude from (I2.28H that the imaginary part of V has 
fixed sign on the interval (oo, 0]. Thus Corollary 12 . 51 shows that y(u) is bounded on the 
interval (— oo,0]. Hence the corresponding Sturm-Liouville solution (j>(u) = exp( f y) 
has no zeros on this interval. In particular, we conclude that ^ (f)(0) = 2t(^o)- D 

Remark. Plotting the function 2l(z) with Mathematica, one sees that on the ray 
e~ J7r / 3 ]R + , the function 21 does indeed have many zeros (and probably an infinite 
number). 



3. Estimates for Parabolic Cylinder Functions 
We now consider a quadratic potential 

V(u) = a + j (u-j) 2 with q,/3, 7 gC. (3.1) 

The corresponding differential equation (jl.ip can be solved explicitly in terms of the 
parabolic cylinder function, as we now recall. The parabolic cylinder function, which 
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we denote by U a (z), is a solution of the differential equation 

*C(*) = (j +a)Va{*). (3-2) 

Setting 

OL 1 

(p(u) = U a {z) with a = — , z = fi* (n-j) , 

a short calculation shows that <j) indeed satisfies (jl.ip . 

According to [8, § 12.5(H)], particular solutions of (|3.2|) can be written as contour 
integrals, 

U a (z) = e~T y e 2 *"# tfH . 

We consider the special solution obtained by choosing T as the contour R + i, where 
we take the usual convention i a_ 2 = exp((a — ^)logt), and the logarithm has its 
branch cut along the negative real axis. It is verified by a direct computation that U a 
satisfies (13. 2ft ; namely, 




where in the last line we integrated by parts. 



3.1. WKB Estimates. Integral representations for parabolic cylinder function have 
been used to derive asymptotic expansions (see for example [7] and [9]). Our goal is to 
derive rigorous estimates with quantitative error bounds which go beyond [7J |9] and 
show that in a certain parameter range, the function U a {z) is well-approximated by a 
WKB wave function. More precisely, the WKB wave function (jl.2p becomes 

<£wkb(*) = ( z 2 ^ 4q ) exp Q/ VC 2 + 4a d(^J 

(z + V^ 2 -6 + 2)5-3 — (3 ' 3) 

= r exp - V ^ 2 - 6 + 2 , 

( z 2_ 6 + 2 )^ F U / 



where we introduced the abbreviation 



1 

a ~2 



We write the parabolic cylinder function as 

U a (z) = e-4 / e /W di , (3.4) 



r 
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where the function / and its derivatives are given by 

f(t) = zt-j-~\ogt, f( t ) = z-t-j t (3.5) 

The zeros of /' are computed to be 

t± = ^(z± Vz 2 - ft) • 

We choose to equal to either t + or t_ such that the angle between z and V z 2 — b is at 
most ^, i.e. 

Re (±z \/z 2 - ft) > . (3.7) 

Then 

± Vz 2 -b 2 = \z\ 2 + \z 2 - b\ + 2Re (±zV z 2 - b) > \z\ 2 + \z 2 - b\ (3.8) 



We now use a "stationary phase-type" approximation to U a in order to derive error 
bounds. To this end, we let / be the quadratic Taylor approximation of /, 

f(t) = /(to) + ^ (t - to) 2 • 
It is convenient to introduce the parametrization 

t = t (l + r); (3.9) 

then 

f(t) - /(to) = -f r 2 + b - (r - log(l + r)) (3.10) 

/»-/(M=(-| + ^)r 2 . (3.11) 
The coefficient of r 2 in (|3.10p is given by 

c--=- t j = -l(z±Vz^b) 2 . (3.12) 

Using (13. 8j) . we obtain 

\c\>l(\z\ 2 + \z 2 -b\) . (3.13) 
Deforming the contour T, we choose the contour to be the straight line 

T {u) = e ia u with u € R and a G [0, vr) . (3.14) 

Then 

:= r(n) - log (l + r(«)) = e ta u - log(l + e ia u) 

4-(e ia s-log(l + e JQ s))(is= / ( e JQ — ) ds 

ds\ &y ') J V l + e ia s) 

u e 2ia s ru e 2ia g n + e -ia g \ ru & 2ia g + e ia g 2 

ds = : 77; ds = - k — ds 



l + e ia s Jo \l + e ia s\ 2 Jo l + s 2 + 2scosa 

e 2ia h{u) + e ia l2(u) , 
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where I\ and I 2 are the real-valued functions 

v ; J l + s 2 + 2scosa 
Choosing a polar decomposition of b and c, 



ds 



1,2 



(3.15) 



6 

8^ 



6=|6|e^, c=\c\e^ with /?, 7 G [— vr, 7r) 
and introducing the abbreviation 

we obtain 

Re (/(*)- /(«o)) = |c| , 

where 

h(u) := cos(2a + 7) u 2 + 2p ( cos(2a + B) Ji(u) + cos(a + B) I 2 {u) 
Lemma 3.1. The function h satisfies the inequality 

P 



(3.16) 



(3.17) 



h(u) , . sin(a + 8) 

-V " cos 2 « + 7) + P — 7T-. — 

2 sm a 



7t z 



< 



2 sin a 



Proof. It is obvious from ([337) and (|3T5jl that /i(0) = = h'(0), and that h"(0) exists 
and is finite. Suppose that u > 0. Then 



/i'(u) < it sup 

s<= 



and integrating from to u yields 

h(u) 



< sup 

seM\{o} 



h'(s) 



h'(s) 
2s 



Similarly, one can estimate h(u)/u 2 from below by inf h'(s)/(2s). Proceeding analo- 
gously in the case u < gives 

p fc'(a) fc(u) h'(s) 
inf < -V < sup 



sSK\{0} 2s u se 

By differentiating (|3.15p . we readily obtain 



2s 



h'(s) ,„ . cos(2a + 5) + s cos(ct + /3) 

— cos(2a + 7) + p - 



2s 



1 + s 2 + 2s cos a 



(3.18) 



which implies that 



h'(s) 

hm — — = cos (2a + 7) . 



s— >±oo 2s 

As a consequence, the function h'(s)/(2s) is bounded from above and below by its 
inner maximum and minimum, respectively. Differentiating (|3.18p and computing the 
zeros, a straightforward computation gives the result. □ 
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Lemma 3.2. By appropriately choosing a in (|3.14p . one can arrange that 



sin a > 



and 
where 



Re(f(t)-f(t ))<-K\c\ |r| 2 , 



K = < 



sin 



1-p 



2 P 



7 

COS ^ 



, 2vr 2vr\ 



(3.19) 
(3.20) 

(3.21) 



cos ( 7 ± 



7T 



2p 1 - sin [p ± 



if 7 



2?r 2tt\ 

i'ii 



(and one can choose the cases ± as desired). If p < 1/8, we can choose k = 1/4. 

3 ' 3 



Proof. In the case 7 £ (—%,%), we choose 



o 



7T — 7 



which gives cos (2a — 7) = —1. Applying Lemma 13. II and using the identities 



sin 



7T — 7 



2 1 a + (3 vr 



1 — sin(a + fi) = 2 sin 



gives (f3T20D . 

In the second case, we apply Lemma 13.11 choosing 



2 sin 2 



P 7 



7T 

a = — 
6 



or a 



5-7T 



A short calculation again gives (|3,20p . 

If p < 1/8, another short calculation shows that k as given by (|3.2ip is always larger 
than 1/4. Therefore, we can choose k = 1/4. □ 

Theorem 3.3. Suppose that the phases f3 and 7 are such that the parameter k in 
Lemma \3.2\ is positive. Then the parabolic cylinder function can be approximated by 
the function 



U a ( z ) = e - £ r+ i «+/(*o) t 



V-8c-6 



with the relative error bounded by 



U a (z) 



U a (z) 



< 



V\8c + b\\b\ 

K 2 |c| 2 



1 + log 2 



( k \c\Y- 

\b\ 



(3.22) 



(3.23) 



Note that the right hand side of (|3,23p tends to zero as \z\ — > 00 in view of (|3. 13[) . 
Proof of Theorem\EM From ([3T4"]) . (^9|) and ([BTLl]) we obtain 

2 poo 

U a (z) = e ~^ +ia+f{to) t / e /(*)-/(*") du . (3.24) 



Replacing / by /, we obtain a Gaussian integral, which can easily be computed to 
obtain ([322]). 
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In order to estimate the integral for large \u\, we fix a parameter L > 0. Then 



/ 



J(t)-f(t ) du 



-L.L] 



e Re(/(t)-/(t )) du ^ 



-L,L] 



-«|c|« 2 du 



< e -5 K l c l" 2 



e 2 



2vr 



e 2 



| K |c| i 2 



K C 



Since / is the quadratic approximation to / at r = 0, the inequality (|3.20p also holds 
for /. Hence the last estimate is also true for /, 



/ 



J(t)-f(to) du 



-L.L] 



< 



2tt 



K C 



e -\n\c\Ll 



Next, on the interval [— L, L] we use the estimates 

- 1 d 



e f(t)-f(to) _ e f(t)-f(to) 



ds 



S f(t)+(i-s)f(t)-f(t ) dg 



(/(*)-/)(*) e -/(*)+d-)/(*)-/(*o) ds 
< ifU\ _ fay e Re (sf(t)+(l-s)f(t)-f(to)) 

<\f{t)-f{t)\e-*M*<\f(t)-f(t)\. 

Moreover, as / is the quadratic approximation to /, we know that 

l-rl 3 , s JXTot | T |3 9 

|/(t)-/(t)|<^su P |/ (3) (r)| < LLn sup ___. 

Since the distance from our contour (|3.14p to the point —1 is equal to sin a, we can 
apply (|3.19p to obtain 

\f(t)-f(t)\<hb\\r\ 3 . 



We conclude that 

i-L 



-L 



;/(*)-/(*0) _ e f(t)-f(t )\ du 



<-\b\ / M 3 du = - \b\L A 



(3.25) 



Combining the above estimates, we conclude that for any L > 0, the following 
inequality holds, 



J(t)~f(t ) _ e f(t)-f(t )\ du 



< 2 ^L e -h^\t 2 + i\b\L*. 



K C 



(3.26) 



We want to choose L in such a way that the two error terms are of comparable size. 
To this end, we set 



L 



2 i i 

lOg 2 



K C 



(K]cj) 

\b\ 



(3.27) 
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< 2 



2tt \b\ 1 ... 4 

At C )2 J K l c l 



K C 



log 2 



.« c a 



|6| 



k 2 |c| 2 



(3.28) 



Computing the numerical constants and using (|3.24|) and (I3.22|) gives the result. □ 
Theorem 3.4. Under the assumptions of Theorem \3~. 



K{z) 



U a (z) 



< 



\b\ 2 W(z)\ 



k 4 |c|V|8c-6| 



1 + log 2 



(K\ C \y 



\b\ 



+ ^^V^^[l + log 3 



lfi|3 



Proof. Introducing the abbreviations 



we have 



/(*) 



d ^„(z)\ d I'(z)-f(z) V{z)_ I{z) - I(z) 

I{z) I(z) I{z) 



dz \U a {z)J dz \I(z), 
Computing the Gaussian integral gives 



8vr 



-8c- b J 

and a straightforward calculation using (|3.28p yields 



?{z) i{z) - I(z) 



I(z) I(z) 
Moreover, 



< 



V2 |6p 



3^ |c| 4 ^|8c + 6| 



\c'(z)\ 1+loe 



(ftjcj) 

l&l 



- ?(*) = / c'(z)r( U ) 2 fe/(*W(*o) - e /(*)-/(*o) , j du 



and thus 



|/'(z)-/»|<|c'(z)| 



2 f e /w-/(* ) _ e m-Ht )\ du 



Proceeding as in the proof of Theorem 13. 3} we obtain similar to (I3.26|) 

u 2 ( e m-m) _ e m-m)\ du 



< 5 e 2 K|6|L +-|6|L 6 . 



(«|c|)3 



Again choosing L according to (|3.27p . we obtain 

1*1 



u 2 Lm-nto) _ e m-f( t0 )\ du 



< 



k 3 |c| 3 



2V2^ + — log 3 



\b\ 



Computing the numerical constants and combining all the terms gives the result. □ 
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We finally explain the connection of the parabolic cylinder function to the WKB 
wave function (13.31): 



Remark 3.5. Assume that \z\ 2 > 4|6|. Choosing to = t+, we know from f|3. 12|) that 



Vz 2 - b 



< \b\ sup d c (y/z 2 -C) < \b\ sup - 

4|C|<|«| 2 4|d<k| 2 2 



Ic| > - \2z 



2 \zV 

> — 
- 2 



2z 



> 



EES) 161 1 

' < 2^8' 

and moreover (|3.7p holds. Hence choosing k = 1/4, Theorems 13.31 and 13.41 apply, show- 
ing that U a (z) is well-approximated by U a {z) with a small error. Next, a straightfor- 
ward calculation yields 



U a = 23 e I +ia v^F 



(z + Vz 2 - b) 
(z 2 -6)1 



1 _ 6 

2 4 



exp I - a/z^ 



(3.29) 



Thus we see that replacing z 2 — b by z 2 — b + 2, we obtain up to a constant, the WKB 
wave function f|3.3|) . Assuming furthermore that |z| is sufficiently large, we conclude 
that U a (z) is indeed well-approximated by 4>wkb(z)- 

In particular, this last observation together with ()3.29j) enables us to compute the 
following limit. 



Corollary 3.6. For any zq £ C, 



lim 



1 U^tZp) = £o 
£+oa t U a (tz ) 2 

3.2. Ruling out Zeros of the Parabolic Cylinder Functions. In analogy to Sec- 
tion [2]3l we now analyze the solution U a {z) along the straight line (|2.25|) . Thus set- 
ting <p(u) = U a (z(u)), this function is a solution of the Sturm-Liouville equation f 1 1 . 1 j) 
with 



+ a 



V{u) = A 2 (i£o+M! 

After suitably choosing A, one can apply Corollary 12.51 to obtain the following result. 

Theorem 3.7. The function U a has no zeros in the complex plane. 

Proof. We always choose A such that A 4 is real. We then obtain 

1 



Im V(«) = - Im [2X 3 z u + A 2 (z 2 + 4a) J . 

Moreover, according to Corollary 13.61 we know that the function y{u) := (p'(u)/4>(u) 
behaves asymptotically for large negative u as 



y(u) = A 



U a (z{u)) 



A 



zq + An 
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We first consider the case KezQ ^ 0. Choosing A = ±i, we get 
Imy(-u) ~ ±~ Rez 

Im V (u) = i Im ( =F 2iz w - (zq + 4a) ) = =p^ Re z - Im + a^ . 

Thus in the limit u — > — oo, the imaginary parts of the functions y and V are non- 
zero and have the same sign. Moreover, by choosing the sign ± appropriately, we can 
arrange that the imaginary part of V does not change sign on the interval (— oo,0]. 
Exactly as in the proof of Theorem 12.61 it follows that U a {zo) ^ 0. 

In the remaining case Rezo = 0, we choose A = a(l + it) with a, r 6 {il}- 
Then A 2 = 2zr and A 4 is indeed real. It follows that 

Imy(u) ~ ^ Im(A 2 ) = ru 
4 Im V(«) = Im ^2A 3 z u + A 2 (zg + 4a) ) 
= 2t Re ^2Az u + (zq + 4a) 
= -4<TnIm(z ) + 2r Re (zq + 4a) . 

Thus in the limit u —> — oo, the imaginary part of y is non-zero, and by choosing a 
appropriately we can arrange that it has the same sign as the imaginary part of V. 
Moreover, by adjusting r we can arrange that ImV does not change sign on the 
interval (— oo,0]. Again proceeding exactly as in the proof of Theorem 12. 6[ we obtain 
that U a (z ) / 0. □ 



3.3. The Airy Limit of Parabolic Cylinder Functions. In this section, we iden- 
tify the parameter range where the parabolic cylinder function is well-approximated by 
the Airy function. The idea is to expand the function / in (|3.4[) around a point to where 
the second derivative vanishes. Then the Taylor polynomial of / around to involves a 
linear and cubic term, just like the exponent in the integral representation (I2.4p of the 
Airy function. The zeros of f"(t) are computed to be ±v6/2. In order to fix the sign, 
we let c be a solution of the equation 




(3.30) 



and set 

t = c 3 . 

Then the function / has the Taylor expansion 

f{t) = /(t ) + [z- 2c 3 ) (t - t ) - -L (t - to) 3 + 0((t - to) 4 )- 

Defining / as the corresponding Taylor polynomial and introducing the parametriza- 
tion 



t = to — cr , 



(3.31) 
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Figure 3. Admissible regions for the integration contour. 



we obtain 

/(*) " /(*o 

fit) ~ f(to 
Then writing (j3.4j) as 



) = (c 4 — cz) T 



f2c 4 



r 2 

-T 2 

2 

r 3 



c b log 1 



cz)t + — 



U a (z) 



( e /(*)-e/W)dt, 



(3.32) 



e /(t) (it + e" 

'r 

the first integral is just of the form of the Airy function (|2.4|) . whereas the second 
integral is the error term. In order to show that the parabolic cylinder function goes 
over to the Airy function, we must make sure that the contour T is compatible with 
that in (12. 4j) . and that the functions e^*** and both decay exponentially on both 
ends of the contour. 

We now demonstrate how to continuously deform the original integration contour t G 
R + i in such a way that the real parts of both / and / are both negative at both ends 
of the new contour. In our parametrization (I3.3ip . we need to specify the contour in 
the variable r. As only the product cr enters in (13.311) . we can assume without loss of 
generality that 

3vr 137r 

— < arg c < mod Zir . 

4 - ° 12 

Then we can deform the contour such that in Figure[3]the integration begins in region B 
and ends in region C (we could also end in region A, but our choice fits well with the 
contour chosen for the Airy function in (|2.4p ). For the real part of — c 2 r 2 to be negative 
at both ends of the contour, we need to assume that (see Figure E|) 



57T 7T 

< 2 arg c < — 

12 6 12 



mod 27T . 



Provided that this inequality holds, there is indeed a contour such that the real parts 
of both / and / are negative at both ends of the contour (see Figure [3|) . Then the 
error term in ()3.32j) decays exponentially on both ends of the contour. 
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Having accomplished this exponential decay, rigorous error estimates can be ob- 
tained in a straightforward way using methods similar to those in Sections 12.11 and 13.11 
We leave the details to the reader. 
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